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Comment on "T-dependence of the magnetic pen- 
etration depth in unconventional superconduc- 
tors at low temperatures: Can it be linear?" Re- 
cently, Schopohl and Dolgov (SD) raised the possi- 
bility that a d-wave superconductor might be thermo- 
dynamically unstable in the limit T — )■ 0. Within the 
framework of linear response theory, they calculated the 
entropy of a d-wave superconductor due to the presence 
of a current, 
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where c/47rA2(q, T) is the (linear) static current response 
function, and B the magnetic field in the sample. They 
argued that the Nernst theorem limT^o S" — > would 
then require dX{q,T)/dT — 0, whereas the penetration 
depth A(r) of a d-wave superconductor is known to be 
linear in temperature, apparently yielding 5(r ^ 0) = 
const. While there are several effects which can lead to 
a nonlinear T-dependence in a superconductor with line 
nodes at low temperatures, including impurities and 
nonlocal effects, [|[^ SD pose an intriguing question of 
principle: Suppose one considers a pure 2D system in a 
special geometry where the nonlocal effects such as that 
discussed in Ref. |Q vanish. In such a case it appears as 
if (5A ~ r. A constant entropy at T — > would therefore 
result, implying the breakdown of BCS theory for a pure 
d-wave state, an unsettling prospect. 

One answer to the question posed by SD has been pro- 
vided by Volovik, [|j who points out that the linear re- 
sponse calculation of SD is not complete, since any phys- 
ical experiment will be performed with a finite current. 
He argues that therefore, when considering Nernst's the- 
orem, the limit T — > should be taken first. Based on 
a semiclassical analysis of the resulting nonlinear effects, 
he shows that the free energy varies as const + rather 
than T at low temperatures, leading indeed to zero resid- 
ual entropy as required for stability. 

While we agree with Volovik's point that nonlinear ef- 
fects always prevent an instability, we believe BCS theory 
is even more robust than implied by his argument, and 
that the answer to the stability question should not rely 
sensitively on the order of the limits taken. As one lowers 
the temperature, there are several ways in which the d- 
wave system can escape the SD "paradox" , of which the 
Volovik scenario is only one. A more general analysis for 
a given particular system with or without a surface should 
include consideration of both nonlinear and nonlocal ef- 
fects, as performed recently for the vortex state. Thus 
in a pure d-wave superconductor, there are three relevant 
energy scales: T, E, 


nonioc 


and E, 


The magnitude 
of the last two energies depend on the particular spa- 
tial variations of the currents in the system (they are not 
independent in the vortex lattice case). In a Meissner ge- 
ometry with generic surface, Enonioc — vj? - q ~ AqS^q/Xq, 
where Aq is the gap maximum, and and Aq are 
the zero-temperature coherence length and penetration 


depth, respectively. Enonioc will be larger than the typ- 
ical nonlinear scale Enoniin — '^f, where is the 
superfiuid velocity, if the applied current is sufficiently 
small. In such a case, the response function entering (1) 
may easily be shown to be X~^{q,T) ^ const + for 
T <^ qvp (only after integrating over all Fourier com- 
ponents does one obtain the asymptotic nonlocal result 
6X ~ r2 of Ref. 0). Thus the "paradox" is trivially 
avoided even in linear response EnonUn 0. 

As discussed by SD, however, in a special situation 
with H _L c II 71 in a half-space geometry, it would seem at 
first that there are no nonlocal effects, since for a strictly 
2D order parameter confined to the ah plane, q • vp = 0. 
In this case, however, the magnetic field is not screened 
(A — > oo) 101 . For a quasi-2D system, on the other hand, 
(\ - wp — qvp ^ 0, and thus nonlocal effects always exist. 
The difference Ck+q/2 — Ck-q/2 (Ck = fc^/2TO — /i is the 
quasiparticle energy in the normal state) which cuts off 
the nodal singularity in the full response also contains a 
further term j^m such that in the SD special case the 
energy scale Enonioc - max[AoCoc/Ao, l/(2mAg)], where 
^Oc is the c-axis coherence length. 

Therefore there appears to be no geometry in which 
nature does not find a way, at sufficiently low tempera- 
ture, to force the penetration depth, and hence the free 
energy, away from a linear T variation, thus avoiding the 
SD "paradox" . A d-wave system may be stable down to 
T = 0. 

P.J. Hirschfeld"'^ M. -R. Li^ and R Wolfle' 
" Dept. of Pliysics, U. Florida, Gainesville, PL 32611 USA 
^ Institut fiir Theorie der Kondensierten Materie, Univer- 
sitat Karlsruhe, 76128 Karlsruhe, Germany 

Received June 1997 

M.-R. Li and P.J. Hirschfeld are grateful for support from 
the A.v. Humboldt Foundation. 


[1] N. Schopohl and O. V. Dolgov, Phys. Rev. Lett. 80, 

4761(1998). 
[2] F. Gross et al, Z. Phys. B64, 175 (1986). 
[3] W.O. Putikka, P.J. Hirschfeld, and P. Wolfle, Phys. Rev. 

B 41, 7285 (1990). 
[4] I. Kosztin and A.J. Leggett, Phys. Rev. Lett. 79, 135 

(1997). 

[5] G.E. Volovik, cond/mat 9805159 

[6] M.H. S. Amin, I. Affleck, and M. Franz, 

cond/mat/9712218 
[7] D. I. Glazman and A. E. Koshelev, Sov. Phys. JETP 70, 

774(1990). 


1 


